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ABSTRACT 

We have investigated the FOV star V398 Aurigae (= 9 Aur) under the assumption 
that it is undergoing non-radial gravity mode oscillations and that the two principal 
periods given by Krisciunas et al. (1995) are correct. We find that the two periods are 
manifestations of an = 3, |to| — 1 spheroidal mode and its toroidal corrections due 
to the rotation of the star. As far as we know, this is the first detection of toroidal 
correction terms in a real star. The two modes probably are the result of rotational 
splitting. 

Our analysis provides for the first time a physical explanation of certain charac- 
teristics of the observed behavior of the star. The amplitude of the radial part of the 
pulsation for /i = 0.795 d~^ is a factor of 4 larger than the one for /2 = 0.346 d~^. 
Since the photometric variability is determined mostly by temperature variations, 
which in turn are determined by the radial part of the pulsation, the photometric 
variability is dominated by the mode with frequency /i. On the other hand, /2 is 
the more pronounced one in all three spectroscopic moment variations (including the 
radial velocity), refiecting that the transverse displacement of /2, and not the one of 
/i, dominates the velocity behavior. 
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Krisciunas & Handler (1995) have compiled a list of 17 stars 
of similar spectral type and luminosity class which appear to 
constitute a new class of pulsating variable stars. These stars 
are typically of spectral type FO to F2, and they are found 
on, or just above, the main sequence in the Hertzsprung- 
Russell Diagram. The best studied examples are 7 Doradus 
(Balona, Krisciunas & Cousins 1994), 9 Aurigae (Krisciu- 
nas et al. 1995), HD 224638 and HD 224945 (Mantegazza, 
Poretti & Zerbi 1994). Given that 7 Doradus is the brightest 
member of the list, and that it was the first one found to 
be variable (Cousins & Warren 1963) , it has been suggested 
that these stars be known as "7 Doradus stars" . 

The photometric time scales of the 7 Dor stars are 0.5 
to 3.5 d - variations an order of magnitude slower than the 
fundamental radial pulsation period for stars of this density. 
9 Aur and 7 Dor have shown evidence for radial velocity and 
spectroscopic line profile variations. In the case of 9 Aur - 
recently designated V398 Aur in the latest named list of new 
variable stars (Kazovarets & Samus 1995) - the radial ve- 
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locity variations appear to follow only one of the two pho- 
tometric frequencies (/2 = 0.346 d~^) found by Krisciunas 
et al. (1995). 

Given the time scale of variations, the evidence for ra- 
dial velocity and line profile variations, and the lack of other 
viable explanations (including star spots), the evidence is 
strong that 7 Dor stars are exhibiting non-radial gravity 
mode pulsations. This has been a surprising observational 
development, because until now only much hotter stars have 
shown evidence for non-radial g-modes (Waelkens 1991). 

Because the power spectrum of the radial velocities of 
V398 Aur did not show evidence for the highest peak in the 
power spectrum of the photometry (/i = 0.795 d~^), Krisci- 
unas et al. (1995) suggest that /2 arises from an ^ = 1 or 
2 spherical harmonic, while /i arises from a higher degree 
harmonic. In this paper we apply the moment method de- 
scribed by Aerts et al. (1992) to the line profiles of V398 Aur 
in an attempt to identify its pulsational mode(s). 



2 CROSS-CORRELATION PROFILES 

Krisciunas et al. (1995) present 95 cross-correlation profiles 
obtained by R. F. Griffin with the Haute Provence Coravel 
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(Baranne, Mayor & Poncet 1979). While the Coravel is used 
primarily for the determination of the radial velocities of 
stars, the cross-correlation profiles give us, effectively, a mea- 
sure of the mean spectroscopic line depths and a parameter 
characterizing the width of the lines. The width parame- 
ter can be interpreted in terms of the projected rotational 
velocity v sin i of the star (Benz & Mayor 1981). The cross- 
correlation profiles can be thought of as the mean line pro- 
files of the star being studied. 

The Coravel data reduction package also gives us a 
Gaussian fit to the points of the cross-correlation profiles. 
In order to carry out the identification of the pulsational 
mode(s) of a star, it is necessary to have a very high sig- 
nal to noise ratio in spectroscopic parameter space. We 
found it necessary to use the mathematical fits to the cross- 
correlation points for our analysis here. (See the sets of 
curved lines in Fig 10 of Krisciunas et al. 1995.) 

3 DETERMINATION OF PULSATIONAL 
MODE(S) 

Mode identification is currently often obtained by means of 
spectroscopic analyses. The moment method used here was 
first introduced by Balona (1986) and was subsequently gen- 
eralized by Aerts et al. (1992) and by Mathias et al. (1994) 
in the case of respectively a mono- and a multiperiodic pul- 
sation. It is based on the time variations of the first three 
moments of a line profile. The periodograms of the three 
moments can immediately be interpreted in terms of the 
periods and amplitudes of the non-radial pulsation (NRP) 
parameters. The observed moment variations are compared 
with theoretically calculated expressions for these variations 
in the case of various pulsation modes. Mode identification 
is obtained by means of a so-called discriminant, which cap- 
tures the discrepancy between the observed and the calcu- 
lated amplitudes. 

The discriminant used in this paper is a generalization 
of the one proposed by Aerts et al. (1992) in the sense that it 
takes into account all seven observed amplitudes of the three 
moments (instead of only three) and it gives them a weight 
according to their uncertainty. This discriminant turns out 
to be more accurate than the one presented by Aerts et al., 
especially in the case of non-axisymmetric modes (Aerts, 
in preparation). Identification is achieved by means of the 
minima of the discriminant, which, in this paper, are defined 
by: 

7^ = min r™(t^p, i, a), (1) 

Vp,i.a 

where Up is proportional to the pulsation velocity amplitude 
(the proportionality constant is given in Aerts et al. 1992), 
i is the inclination angle, and a is the width of the intrinsic 
profile that is assumed to be Gaussian. The best solution 
in the space of the pulsation parameters {£,m,Vp,i,a) is 
the one with the lowest minimum 7™. We recall that the 
discriminant is unable to find the sign of the azimuthal 
number m, and that it is an accurate identification tool 
for low-degree modes {£ < 4). In constructing the discrimi- 
nant given by equation we have assumed that the v sin i- 
value of V398Aur is accurately known and adopt 17.8 km 
s~^ (Krisciunas et al. 1995). However, since these authors 
claim that the non-radial pulsations of the star are giving 



rise to line width variations, they point out that the "true" 
projected rotational velocity may actually be smaller than 
the mean observed mean value of 17.8 km s~^, perhaps as 
small as the smallest observed values of the line width (~ 
15 km s"^). 

We here apply the biperiodic moment method to the 
profiles of V398 Aur. The three observed moment variations 
are shown as dots in Fig|l|, where they are put in phase 
according to /i (right) and /2 (left) respectively. The peak- 
to-peak amplitude of the radial velocity (=first moment) 
amounts to 3.0km s~^ for /i and to 5.2km s~^ for /2. It 
is clear that a frequency analysis of the radial-velocity data 
cannot yield the same accuracy as the photometry presented 
in Krisciunas et al. (1995). Nevertheless, we performed a pe- 
riod analysis on the three moments of V398 Aur to check if 
/i and /2 are present in them. We used the Phase Dispersion 
Minimization method of Stellingwerf (1978) and found that 
the frequency /2 is present in all three moment variations, 
especially in the odd ones. The presence of the frequency /i 
is less clear, but beat- and sum- frequencies of /i and /2 are 
present in the second and third moment. We further studied 
the variations of the moments assuming that the two fre- 
quencies /i = 0.795 d~^ and /2 = 0.346 d"'^ are accurate. 
The solid lines shown in Fig|l] are the theoretical fits to the 
observations for a model with /i (right) and /2 (left). The 
fact that the third moment is centered around km s~^ re- 
flects that there is no variability on a time scale larger than 
the time span of our data. 

We can determine the ratio of the horizontal to the ver- 
tical pulsation velocity amplitude K by means of the bound- 
ary conditions as a function of stellar mass and radius if 
we know the pulsation frequency in the corotating frame. 
Since we only know the frequencies in the observer's inertial 
frame, we use these to approximate the true K-vahies. The 
actual mass of 9 Aur A, the primary, is not known. While 
the B component is an M2 V companion at a distance of 
~100 an (Krisciunas et al. 1993), the orbit of that compan- 
ion is not well known enough to determine the mass func- 
tion of the system. Adopting a mass range of 1.5-1.7 Mq, 
which is reasonable for FOV stars (Popper 1980) , and further 
R = 1.647?0 (Mantegazza et al. 1994), we find Ki G [40,45] 
for /i and K2 € [211,240] for /2. It is clear from these 
numbers that we are dealing with high-order g-modes: the 
transverse displacement is dominant over the radial one. 

The equivalent width of the profiles varies by 8% for 
/i and by 3% for /2. Since equivalent width variations are 
usually interpreted in terms of temperature variations, these 
findings are in agreement with the photometric variations in 
the sense that the most important temperature variation is 
due to /i. 

It is clear from Fig|l] that some of the theoretical am- 
plitudes of the moments are uncertain. For such data the 
discriminant defined in equation (|l|) is much better than the 
one defined by Aerts et al. (1992), because the former is con- 
structed in such a way that the most accurate amplitudes 
are given the largest weight. The amplitudes of the terms of 
the second moment varying with frequencies /i and /2 are 
at least as large as those varying with frequencies 2/i and 
2/2 . This is a signature of the presence of non-axisymmetric 
modes (see Aerts et al. 1992). 

The amplitudes of the moments are used to calculate 
the discriminants r™(up, i, a) and their minima 7™ for each 
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Figure 1. The first, second, and third moment of the mathematical fits to the cross-correlation profiles of V398 Aur taken between 1993 
December 25 and 1994 January 10 UT (see also in Section|2|). The dots represent the observations, while the full lines are the fits for a 
model with /2 (left) and /i (right). 



Table 1. The different minima of the discriminants for the two modes of V398 Aur. 7™, tip, and cr are given in km s^^. The radial 
component of the velocity has an amplitude proportional to Dp, while the transverse component's amplitude is proportional to Kvp (see 
also Discussion). 
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candidate mode {£,m). We have used the upper hmits for 
the if- values as input numbers for the discriminant. In the 
case of V398 Aur, the solution is independent of this choice 
since the discriminant searches the most likely amplitude 
of the transverse velocity Kvp-. a lower input K gives a 
higher Wp-value and vice versa. The adopted K-valne is im- 
portant when the radial and transverse velocities are of the 
same order of magnitude. The results of the discriminant 
for vsini = 17.8km s~^ are listed in Table^ We only give 
the best solutions in parameter space. It is seen from the 
table that the two frequencies give rise to an almost iden- 
tical list of most likely modes. This situation is comparable 
to the results found for the /3 Cephei star P CMa (Aerts et 
al. 1994) and are interpreted in terms of the presence of 
two modes with identical degree £ but opposite azimuthal 
numbers m. The combination of two £ = 3 modes, one with 
m = +1 and one with m = —1, seems most likely here, 
since the inclination and intrinsic width fits nicely for such 
a combination, while this is less the case for the £ = 4 or 2 
solutions. Nevertheless, we cannot exclude the possibility of 
other combinations of modes given in Table |l], since the er- 
ror on the inclination can be quite large. It does seem fair to 
conclude that |m| — 1 for both modes. We have recalculated 
the discriminant for the lower limit of the projected rotation 
velocity, i.e. usini = 15 km s"'^. The results are comparable 
to the ones given in Table |l[ but with slightly different values 
of the inclination, the amplitude, and the intrinsic width. 

The question then arises if the two modes can be due to 
rotational splitting. Dziembowski & Goode (1992, equations 
22 and 117) give an expression for the observed frequencies 
/ expected from rotational splitting of a mode {£, m) in the 
case of rigid rotation and in the limit of high-order (;-modes : 

' - ■^''-"C"^) (2) 
n?^^ U{£ + 1) (2^(f + 1) - 3) - 9 
/o 2£2(^-Kl)2(4£(£+l) -3) ' 

where /o is the frequency in the corotating frame in the case 
of the absence of rotation. This formula is derived under the 
assumption that n « /o- Uniform spacing occurs if the 
second-order term in i7 is negligible. If we assume that this 
is the case, then /o = 0.571 d~^. If the £ = 3,m = +1 solu- 
tion belongs to /2 and the £ = 3,m = —1 one to /i, then we 
obtain this /o from /2 and from /i for Q, = 0.245 d~^. As- 
suming again that vsini — 17.8km s~^ and a stellar radius 
of 1.64/?0 (Mantegazza et al. 1994), such a rotation fre- 
quency corresponds to an inclination angle of 61°, a value 
fully in agreement with the ones derived from the two dis- 
criminants. For fo = 0.571 d~^, the second-order term in 
equation (^) amounts to only 0.008 d~^. None of the other 
combinations of = £2,mi — — m2 listed in Table |l] have an 
inclination that is close to 61° for both modes. We conclude 
that the two modes of V398Aur can be due to rotational 
splitting of an ^ = 3 mode. Unfortunately, the suggestion of 
Krisciunas et al. (1995) that the "true" projected rotational 
velocity may be less than the observed mean value cannot 
be confirmed by a comparison of the results of the discrimi- 
nants for different v sin i- values, because the uncertainty on 
the inclination is unknown. The true ratios of the horizontal 
to vertical velocity amplitudes for uniform splitting of the 
£ = 3 mode are Ki e [83, 94], 7^2 G [72,82]. There is no 



evidence of the presence of other frequency components of 
the uniform rotational splitting belonging to ^ = 3 in the 
photometry. 

Assuming a rotation frequency Q of 0.245 d~^ and a 
frequency fo — 0.571 d~^ as found in the case of uniform 
rotational splitting of the two £ — 3 modes, we obtain 
Q/ fo = 43%. With such a high ratio, the frequency-splitting 
formula given in equation ^ may not be very accurate, but 
it is the best one available. Also, the velocity field of an 
NRP can no longer be described in terms of one {£, m)- 
value, because the Coriolis force and the centrifugal forces 
are not negligible. The Coriolis force induces toroidal cor- 
rection terms described by an and an i7!!i spherical 
harmonic in the case of an NRP with wavenumbers {£, m) 
in the non- rotating case (see e.g. Aerts & Waelkens 1993). 
It is clear that our results found from the moment method, 
which does not take into account these corrections, have to 
be interpreted with caution. It is possible that some of the 
solutions with |m| = 1 found by the discriminant are mani- 
festations of the toroidal corrections since the latter have the 
same azimuthal number as the non-rotation spherical har- 
monic. If this is correct, then it would be, as far as we know, 
the first detection of toroidal corrections in a real star. 

In order to have an idea of the line-profile variations 
(LPVs) expected for the combination of the two £ = 3 modes 
as found by the discriminant, we have generated theoretical 
LPVs with the code presented by Aerts & Waelkens (1993). 
This code takes into account the toroidal correction terms 
due to the Coriolis force by means of a perturbation anal- 
ysis assuming that Sl//o « 1, but not those related to the 
centrifugal forces. As far as we know, no line-profile code 
is available that does take into account the latter effects. 
Because of the large S7//o-value, it is questionable if the 
perturbation analysis is accurate. With this in mind, we do 
not expect perfect fits, the more so since our code does not 
include equivalent width variations and assumes an intrin- 
sic Gaussian that is time-independent. It is clear that these 
conditions are not exactly fulfilled. Furthermore, the effec- 
tive temperature of the hemisphere of the star facing us is 
not constant. From the variations of B-V color of the star 
(Table 3 and Fig 6 of Krisciunas et al. 1995) and the rate 
of change of temperature of mid-main sequence stars as a 
function of B-V (Allen 1973), we estimate for V398 Aur 
that Tcff varies by « 100 K. Nevertheless, we can be confi- 
dent that the solution proposed by the discriminant is good 
if the obtained theoretical LPVs are compatible with the 
observed ones. 

The width of the intrinsic profile found by the discrim- 
inant must be an overestimate of the true value, since the 
velocity contributions of the toroidal terms were neglected 
in the non-rotating model. This width can, however, easily 
be found from the observed LPVs once a mode identification 
has been obtained. We have found that a — 5 km s~^ is a 
more realistic value. Also the pulsation amplitudes will be 
overestimated by the discriminant for the same reason. On 
the other hand, the 7f-values used were also overestimated. 
We have generated theoretical LPVs with all the parame- 
ters except the «p-values fixed (we have again chosen the 
upper limit for both K -values of zeroth-order in SI, but as 
explained above, this choice does not matter). In this way, 
a better estimate of the Up-values could be obtained for the 
two modes. Finally, the phase shift between the two modes 
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for some monoperiodic p-modes but has not yet been stud- 
ied in the case of (;-modes or multiperiodicity. A detailed 
study about this is currently being undertaken (Schrijvers 
et al., in preparation). The fact that the differences between 
the two sets of profiles are not very large makes us confident 
that the results of the discriminant are not too bad. We 
have also calculated theoretical profiles for the other can- 
didate rotationally splitted modes listed in Table|^, i.e. the 
I = A, \m\ = 1 and the 1 = 2, \m\ = 1 solution. In doing so, 
we have taken the average of the two listed i- values. In both 
cases, the correspondence between the observed profiles and 
the theoretical ones is worse compared to the £ = 3, \m\ = 1 
case described by (^). 



Figure 2. A comparison between observed and theoretical LPVs 
for the parameters listed in (^). The dots are actual observed 
points from the cross-correlation profiles obtained by R. F. Griffin. 
The solid lines are profiles based on a theoretical model taking 
toroidal terms into account. The dashed lines are profiles for a 
theoretical model that neglects rotational effects. 

was determined from the first moment. The results of our 
theoretically generated LPVs for the best choice of the pulsa- 
tion amplitudes are shown in Fig^ where we compare them 
with some of the observed profiles. We show 10 arbitrarily 
chosen observed profiles that are spread in phase according 
to /2 with respect to an arbitrary reference epoch. The solid 
lines are the profiles obtained with the code that takes into 
account the toroidal terms {O,/ fo = 0.43), while the dashed 
profiles are those obtained with a model that neglects the 
rotation effect (i.e. those with Q/fo = 0.0). As other input 
parameters, we used : 

/i = 0.795d-S^i = 3,mi = -1, 

Up = 0.100km s-\Ki = 94,Kivl = 9.4km s"^ 
< /2 = 0.346d~\£2 = 3,m2 = +1, (3) 

Up = 0.025 km s~\ K2 = 82, K2vl = 2.1 km s"^ 
vsini = 17.8km s^^,a = 5km s~^,i = 60°. 

v 

First of all, the observed and calculated profiles show the 
same global behavior such that our proposed solution of the 
velocity parameters listed in remains valid. The varia- 
tion in line depth is accounted for by some of the theoretical 
profiles, but not by all of them. The evolution of the two 
sets of theoretical profiles during the cycle are not too dif- 
ferent from each other. It seems that the profiles without the 
toroidal terms are better during the first part of the cycle, 
while the others are better during the second part. The effect 
of the toroidal terms on LPVs depends completely on the 
kind of modes. It was studied by Aerts & Waelkens (1993) 



4 DISCUSSION 

Zerbi et al. (1996) have analyzed the photometric varia- 
tions of V398Aur using data from the 1994/5 observing 
season. While they find essentially the same principal fre- 
quency in the power spectrum of the photometric data (/i 
= 0.796 d~^), they do not find a single peak at /2 = 0.346 
d"'^. Instead they find multiple frequencies in the range 
0.28 < / < 0.49 d~^. It is known that the amplitudes of 
the sinusoids making up the photometric variations are not 
constant (Krisciunas et al. 1995); these authors make a case 
that the frequencies /i and /2 were constant from 1992 Jan- 
uary to 1994 February. The large data set obtained by Zerbi 
et al. in 1994/5 allows us to wonder if the frequencies in 
the power spectrum of the photometry of V398 Aur are also 
variable. In any case, it was found from the spectroscopic 
data of the 1993/4 season that the frequency 0.346 d~^ is 
the most pronounced one in all three moment variations. 

It was emphasized by Krisciunas et al. (1995) that the 
difference in photometric and spectroscopic behavior of the 
frequencies /i and /2 can be understood if /2 would be due 
to a low-degree mode and fi to a higher-degree one. This 
explanation is no longer appropriate if our identification is 
correct, since both modes have the same geometrical repre- 
sentation. From the parameters given in (|^), we note that 
the amplitude of the radial part of the pulsation for /i is 
a factor 4 larger than the one for /2, while its angular de- 
pendence is the same. Since the photometric variability is 
determined most of all by temperature variations and since 
the latter are determined by the radial part of the pulsation, 
it is quite understandable that the photometric variability 
is dominated by the mode with frequency /i . 

For the spectroscopic behavior, the radial displacement 
can completely be neglected in the case of large transverse 
amplitudes as we have here. The relative importance of both 
modes for the line profiles could be the following. The dis- 
placement field ^ correct up to first order in 0./fo in the 
case of the two 1 = 3 modes can easily be calculated, e.g. 
from the paper by Aerts & Waelkens (1993). The toroidal 
corrections for the prograde mode (m = — 1) and for the 
retrograde mode (m = -1-1) are different. Therefore, it is 
possible that the total transverse velocity is larger for the 
mode with frequency /2 than for the mode with frequency 
/i, despite its smaller spheroidal amplitude KA. Quite un- 
derstandably, the line profiles are then mostly determined 
by /2. 

In order to test this point of view, we have explicitly 
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calculated the velocities in the direction of the observer for 
both modes, taking into account the toroidal corrections. 
We indeed find that, although the retrograde mode has a 
four times smaller spheroidal zero-rotation amplitude, its 
total velocity component in the line of sight becomes larger 
due to the toroidal correction terms. The mgiximum velocity 
reached by the retrograde mode is a factor 1.8 larger than 
the one of the prograde mode. 

It thus seems possible that the toroidal corrections that 
appear due to influence of the rotation on a spheroidal mode 
introduce a difi^erent photometric and spectroscopic behav- 
ior when both prograde and retrograde modes are involved. 
To our knowledge, this is the flrst detected consequence of 
the presence of such toroidal corrections. Our point of view 
could be further verified by means of an analysis of a large 
set of high-resolution, high S/N spectra of a star pulsating 
in both prograde and retrograde modes and having a large 
rotation frequency, i.e. having f2//o > 20%. 
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